Traffic and the visual perception of space 
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. Abstract 



During the attempt to line up into a dense traffic people have necessarily to share a limited space under turbulent 
conditions. From the statistical point view it generally leads to a probability distribution of the distances between the 
traffic objects (cars or pedestrians). But the problem is not restricted on humans. It comes up again when we try to 
' describe the statistics of distances between perching birds or moving sheep herd. Our aim is to demonstrate that the 

O . spacing distribution is generic and independent on the nature of the object considered. We show that this fact is based 

^ ■ on the unconscious perception of space that people share with the animals. We give a simple mathematical model of this 

' phenomenon and prove its validity on the real data that include the clearance distribution between: parked cars, perching 

birds, pedestrians, cars moving in a dense traffic and the distances inside a sheep herd. 

^ ' 1 Introduction 

Everyone knows that to park a car in the city center is problematic. The amount of the available places is limited and it has 
to be shared between too many interested parties. Birds face the same problem when a flock tries to perch on an electric 
^ line. Similar reasons lead also to the troubles in the highway traffic, to the pedestrian queues and so on. For instance the 

unpopular transport jams are a consequence of the car-car interaction in a regime of a high car density. 
QQ , The interaction is basically evoked by the brain activity and mediated through the muscles (pedestrians) or accelera- 
' tor/brake pedals (cars). In both cases it is the brain that is responsible for the interaction. So it should be not surprising 
, that the phenomena observed for cars and pedestrians are similar. Even more: the spatial perception is evolutionary very 
old and people share it with animals. So we should be not surprised to find the same results also when dealing with animal 
0^ . herds instead of human traffic. Anyhow: a deeper understanding of the related processes is of interest and will be discussed 
' here. 

J> , The attempt to describe the intrinsic and basically unconscious mechanisms used by the human brain on the basis of 
their everyday activity is not new. It goes back at least to the celebrated work of G.K. Zipf ITI describing the universal 
rS features of languages. The concept was worked out in detail by the sociologist H. A. Simon, within a set of assumptions 
which became known as the SimonSs model |2], see also [3] for the recent results in musicology. 

The subjects of our study (parking, car transport, pedestrian dynamics, herd and flock dynamics) were treated sepa- 
rately in the past. The random car parking model was introduced by Renyi ID (see fS], ||6l for review) and compared with 
the data collected on the street in l?], [S). Other models describe the car transport on highways - see for instance ifTOl and 
||9| for review. Another approach is used for the pedestrian dynamics utilizing social forces - see |[TT|. For the dynamics 
of a starling flock see |[T2l and ITSl 

Our aim is to present a simple theory based on the visual perception that is capable to describe all the observed 
facts regardless whether they originate from cars, pedestrians or animals. The perception mechanism is very old and 
therefore shared by many species ranging from insects to mammals. For humans it is processed automatically and without 
the conscious control. We will use it to understand the statistical properties of the distances between the neighboring 
competitors (cars, pedestrians, birds and sheeps) in a situation when the available space is limited. . 
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The paper is organized as follows: In the Section 2 we describe the mathematical model and derive a one parameter 
family of possible distance distributions. Section 3 contains the psychophysical background of the distance control based 
on the unconsciously evaluated time to collision with the neighbor. The Section 4 contains the comparison of the theory 
with the collected data. 

2 Dividing the available space 

To illustrate the approach we focus on the spacing distribution (bumper to bumper distances) between cars parked in 
parallel. The generalization of the method to other situations will be discussed at the end of this section. 

We will assume that the street segment used for parking starts and ends with some clear and non-transportable part 
unsuitable for parking. It can be a driveway or a turning to the side street. Otherwise the parking segment is free of any 
kind of obstructions. We will assume that it has a length L and is free of any kind of marked parking lots or park meters. 
So the drivers are free to park the car anywhere in the segment provided they find an empty space to do it. We suppose 
also for simplicity that all cars have the same length ^o- Since many cars are cruising for parking there are not free parking 
lots and a car can park only when another parked car leaves. To simplify the further formulation of the problem and to 
avoid troubling with the boundary effects we assume that the street segment under consideration forms a circle. The car 
spacing distribution is obtained as a steady solution of the repeated car parking and car leaving process. 

To park a car of the length one needs (due to the parking maneuver) a lot of a length larger then k, 1.2/o- So in a 
segment of the length L the number of the parked cars equals Xo N k, [L / {1.21q)\. Denoting by Dk the spacing between 
the car k and fc + 1 we get = L — NIq and after a simple rescaling finally 

N 

^i?fc = l. (1) 

k=l 

Since all lots are occupied the number of the parked cars is supposed to be fixed. The repeated car leaving and car 
parking reshuffles however the distances Dk- We will treat them as independent random variables constrained by the 
simplex ([T]). The distance reshuffling goes as follows: In the first step one randomly chosen car leaves the street and the 
two adjoining lots merge into a single one. In the second step a new car parks into this empty space and splits it again into 
two smaller lots. The related equations are simple. If a car leaves the neighboring spacings - say the spacings _D„, D„+i - 
merge into a single lot D: 

D = D„ + an+i+k- (2) 
When a new car parks to D it splits it into £)„, Dn+i'- 

bn = a{D-lo) 
bn+i - {l-a){D-la). (3) 

where a £ (0, 1) is a random variable with a probability density q{a). The distribution q(a) describes the parking pref- 
erence of the driver We assume that all drivers share the same q{a). (The meaning of the variable a is straightforward. 
For a = the car parks immediately in front of the car delimiting the parking lot from the left without leaving any empty 
space. For a = 1/2 it parks exactly to the center of the lot D and for a = 1 it stops exactly behind the car on the right.) 
Combining (|2]i and (O gives finally the distance reshuffling 

Dn = a{Dn+Dn+l) 

bn+i = + (4) 

and the car length Iq drops out. The simplex ([T]i is of course invariant under this transformation. 

For various choices of n the mappings (HI are regarded as statistically independent. Since all cars are equal the joint 
distance probability density P{Di, ...,Dn) has to be exchangeable ( i.e. invariant under the permutation of the variables) 
and invariant with respect to Its marginals pk{Dk) (the probability density of a particular spacing Dk) are identical: 

PkiDk) ^ p{Dk) = [ P{Di,...,DN)dDi..dDk-idDk+i..dDN. (5) 

JDi + .. + Djv = 1 
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A standard approach to deal with the simplex ([TJ is to take as independent random variables normalized by a sum: 

= ^rr- (6) 

where dk are statistically independent and identically distributed. Moreover: it is obvious that the distribution of {Di, ..,Dm} 
is invariant under the transform dU merely when the distribution of {di, .., djv} is invariant. The relation ^ reads for the 
variables (i„: 

dn = a{dn + dn+i) 
dn+1 = [l - a){dn + dn+l) (7) 

(a, dn, dn+1 are now statistically independent and d„, dn+i are identically distributed). 

The parking maneuver is regarded as known and described by the distribution q{a). For simplicity we assume a 
symmetric maneuver, i.e. q{a) = g(l — a). This means that the drivers are not biased to park more closely to a car 
adjacent from the behind or from the front. With given q{a) we look for the distribution of d„ that is invariant under the 
transform (|7|. In other words the effort is to solve the equation 

d^a{d + d') (8) 

where d' is an independent copy of the variable d and the symbol = means that the left and right hand sides of (O have 
identical statistical properties. 

Distributional equations of this type are mathematically well studied - see for instance lfT4ll - although not much is 
known about their exact solutions. In particular it is known that for a given distribution q{a) the equation (O has an unique 
solution that can be obtained numerically by iterations. We are however interested in an explicit result. We restrict there- 
fore the possible densities q{a) to a two parametric class of the standard /? distributions. Then the solution of (O results 
from the following statement ifTSl : 

Statement: Let di,d2 and a be independent random variables with the distributions: di ~ r(ai, 1), d2 ^ r(a2, 1) 
and a ~ /3(ai, a2). Then a{di + ^2) ^ r(ai, !)■ 

The symbol ~ means that the related random variable has the specified probability density. V{g, 1), /3((7i, 172) denote 
the standard gamma and beta distributions respectively. 

For a symmetric parking maneuver gi = 32 = 5- The variables di, ^2 are then equally distributed and the solution of 
([8]) reads d ~ r(g, 1). The relation ^ returns the spacings Dk- We find that the joint probability density P{Di, Dm) 
is nothing but a one parameter family of the multivariate Dirichlet distributions on the simplex ([T]) |[T6l : 



P{D,,...,Dm) = ^,^Dr'Dl-\..D%-' (9) 



Its marginal Q is simply D ~ P{g, {N — l)g). Normalizing the mean of to 1 we are finally left with 

p{D) - 1/3 (g, {N - l)g, ^) . (10) 



A similar reasoning applies also for the moving cars. Assume a sequence of cars in a high density traffic. All the cars 
move with similar velocity and with the mutual distances Dk- In course of the traffic flow the driver of the car k tries to 
optimize his position. He can overtake the car fc + 1 or update the distances to the neighboring cars fc — 1 and fc + 1. In 
doing so he uses the same mechanism as for parking. The only difference is that in the traffic flow we assume the existence 
of certain safety margin. So the update mapping reads now 

Dn = d + a{Dn + Dn+1 - 2d) 
Dn+1 = d+{l-a){Dn + Dn+i-2d). (11) 

where d is the safety margin representing the minimal distance. For maneuvers related with the approach to a standing 
object (like the parking maneuver) we set d = 0. In the course of a car-following in a dense traffic the value of d reflects 
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the reaction time of the driver and his velocity. Changes in the driving situation like night versus day conditions affect this 
component. We will however neglect this fact and regard d as a constant. For experimental results describing the value of 
d under various conditions see for instance ifTTll . 

We assume that outside the safety margin the driving strategy is identical with that of the parking. So in a steady traffic 
flow the distribution of the distances Dk — d has to conforms with the distribution obtained for the car parking. 

The cars are of course not somehow extraordinary. The same reasoning applies for pedestrians, animals in a herd and 
so on. In all cases the distribution q{a) is crucial. We will argue that q{a) is related to the inborn visual perception of the 
distance and independent on the "hardware" actually used to realize the motion. 



3 Distance perception 

The ranging maneuver is described by the probability density q{a) and defined by the relation To ensure a solvability 
we restricted the possible maneuvers to q{a) = /3{g, g, a), with g being a free parameter We will now demonstrate that 
the natural choice is g = 3. 

The point is that for small a the behavior of q{a) reflects the capability of the driver/pedestrian/animal to estimate 
small distances. The collision avoidance during the ranging is guided visually. We assume that the same visual ability 
is shared by all participants. If this applies the behavior of q{a) for small a has to be generic, i.e. independent on the 
particular situation. It is fixed merely by the perception of the distance. 

A distance perception is a complex task and there are several cues for it. Some of them are monocular (linear per- 
spective, monocular movement parallax etc.), others oculomotor (accommodation convergence) and finally binocular (i.e. 
based on the stereopsis). In human all of them work simultaneously and are rehable under different conditions - see ifTsl 
for more details. For the ranging however the crucial information is not the distance itself but the estimated time to collide 
with the neighbor which has to be evaluated using the knowledge of the mutual distance and velocity. 

It has been argued in a seminal paper by Lee |fT9l that the estimated time to collision is psychophy sic ally derived using 
a quantity defined as the inverse of the relative rate of the expansion of a retinal image of the moving object (this rate 
is traditionally denoted as r). Behavioral experiments have indicated that t is indeed controlling actions like contacting 
surfaces by flies, birds and mammals (including humans): see 1201 . 1211 . 1221 . Moreover the studies have provided abundant 
evidence that r is processed by specialized neural mechanisms in the the retina itself and in the brain l23l . The hypothesis 
is that T is the informative variable for the collision free motion - see ||24| for review. 

Let 9 be the instantaneous angular size of the observed object (for instance the front of the car we are backing to 
during the parking maneuver). Then the estimated time to contact is given by 

(12) 



d9/dt 

Since 9{t) = 2 m:ctan{Lo/2D{t)) with Lq being the width of the approached object and D{t) its instantaneous distance, 
we get 

tU) = ^- , arctan f- . (13) 

^' 2Lo{dD{t)/dt) \2D{t)J ^ ' 

For D >> Lq and a constant approach speed v = —dD/dt the quantity r simply equals to the physical arrival time: 
r = D/v. For small distances, however, r w / (vLq) and the estimated time to contact decreases quadratically with 
the distance. (Note that t gives the arrival time without explicitly knowing the mutual velocity, the size of the object and 
its distance.) 

We assume that the probability to exploit small distances is proportional to the estimated time to contact . This means 
in particular that if t evaluated in the course of an approach is small (i.e. a collision is impending) the maneuver is 
stopped. Based on this principle we get for small distances p{D) w and from[3]finally q{a) « a? for small a. Since 
q{a) ~ (3{g, g, a) this sets the parameter g to 5 = 3. The normalized clearance distribution ( fTOl i reads simply 

p{D) = 1/3 f 3, 3(iV .Jjf\.. D\N Or-^ (14) 



iV"^ V \N J 2r(3(7V-l)) 

The described mechanism works so to say in the background, i.e. without being conscious. Moreover: r is evaluated 
equally by humans and by animals. We will show in the next section that this fact leads to an universality in their behavior. 
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Figure 1: The probability density p{D) evaluated for the measured data (crosses) is in a log-log scale compared with the 
function 2.3 * (full line). The agreement for small D is evident. 

4 The measured data 

The estimation of the distance D through t leads simply to p{D) w for small D. So let us first check the validity of 
this relation. There exist a simple observation that enables us to do it: the car stopping on a crossing equipped with traffic 
lights. If the light is red the cars stop and form a queue. We assume that the drivers stop independently and in a distance 
to the preceding car that is evaluated by r. So the clearance statistics should give an evidence of the validity of the r 
hypothesis. There was also a direct experiment measuring the clearance statistics in laboratory conditions - see Il25l . ll26ll . 

We photographed the car queues in a front of the red light. The photographes were taken all from one spot and at the 
same daytime. The clearances were finally obtained by digitalization. Altogether we extracted 1000 car distances from 
one particular crossing in the city of Hradec Kralove (Czech republic) and evaluated the corresponding probability density 
p{D). Similar measurement has been done also on several crossings in Prague - see 1271 . If there is a linear correlation 
between the stopping distance and the estimated time to contact r the obtained distance density p{D) should behave as 
« for small D. The result for small distances is plotted on the Figure [T] It shows a very nice agreement with that 
assumption. 

To show that the r mechanism is generic leading to a distance distribution that is independent on the objects (the 
object can be a man, animal or a car) we divide the further observations into two categories: the first contains sedentary 
objects and the second object moving in a dense environment. 

Let us start with the the clearance distribution obtained for cars parked in parallel and for birds perching on a power 
line. In both cases the "parking segment" is full, i.e. there is not a free space to place an additional participant. We have 
argued that under this conditions the resulting distribution is invariant under the transform ^ with the parameter g in (fTOl l 
fixed to 3. The "parking" segments under consideration were long and containing a large number of objects. So >> 1, 
and the constrain ([TJ does not play a substantial role. In this case p{D) equals to r(3, 1, D). 

To verify the prediction we measured the bumper to bumper distances between cars parked in the center of Hradec 
Kralove (Czech Republic). The street was located in a place with large parking demand and usually without any free 
parking lots. Moreover it was free of any dividing elements, side ways and so on. Altogether we measured 700 spacings 
under this conditions. 
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Figure 2: The probability density of the distances between the parked cars (crosses) and perching starlings (squares) is 
compared with the prediction of the theory (full line). The mean distance is normalized to 1. 

For the birds we photographed flocks of starlings resting on the power line during their flight to the south. The line 
was "full", i.e. other starlings from the flock were forced to use another line to perch. The bird-to-bird distances were 
obtained by a simple digitalization - altogether 1000 bird spacings. After scaling the mean distance to 1 the results were 
plotted on the Figure |2] and compared with the prediction ( fTOl ). 

The probability distributions resulting from these data seems to be (up to the statistical fluctuations) identical and in 
a good agreement with the model prediction. This is amazing since the used "hardware" is fully different. The under- 
lying psychophysical mechanism for the time to contact estimation, is, however, identical. (For the experimental results 
concerning the relevance of r for the space perception of pigeons see |28|.) 

Let us now pass to the traffic streams, i.e. to a situation when the objects are collectively moving. We assume a 
dense traffic. The distances between the moving objects are small and have to be constantly controlled to avoid possible 
collisions. We use the mapping ( fTTI ). 

Outside the safety margin the collision avoidance is assumed to be based on t. This means that the distribution of 
the distances Dk — d follows the mapping (|7| with the same probability density q{a) as in the "parking" case. So the 
distributions of Dk — d should be universal. 

In order to verify this hypothesis we organized a simple experiment with pedestrians in a narrow corridor. Using two 
light gates we measured their velocity and the time interval that elapsed between two subsequent walkers. This enables us 
to reconstruct the mutual distances and evaluate the distance probability density. The same device and method was used 
also for a sheep herd moving through an aisle between two near yards. The third source of data are cars moving on a 
highway in a dense traffic. The velocity and time stamps of the individual cars were obtained by induction loops placed 
below the roadway. 

The safety margins d are clearly different in these three cases and we removed them by subtracting the minimal 
distance from the given data set. The mean distance was finally normalized to 1. The result is plotted on the figure [3] 

Again: the resulting distance distribution is universal and in agreement with the theory based on r hypothesis. 

To summarize we have demonstrated that the clearance between objects (cars,pedestrians,birds and sheeps) is largely 
universal. This surprising observation can be understood as a consequence of an universal distance controlling mechanism 
shared by human and animals. 
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Figure 3: The probability density of the distances between the cars on highway in a high density regime, pedestrians 
walking in a narrow corridor and sheeps moving in an aisle between two barriers is plotted. The minimal distance on each 
particular data set is subtracted and the mean distance is normalized to 1 . 
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